Rotating superfluids in anharmonic traps: 
From vortex lattices to giant vortices 
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We study a superfiuid in a rotating anharmonic trap and explicate a rigorous proof of a transition 
from a vortex lattice to a giant vortex state as the rotation is increased beyond a limiting speed 
determined by the interaction strength. The transition is characterized by the disappearance of the 
vortices from the annulus where the bulk of the superfiuid is concentrated due to centrifugal forces 
while a macroscopic phase circulation remains. The analysis is carried out within two-dimensional 
Gross-Pitaevskii theory at large coupling constant and reveals significant differences between 'soft' 
anharmonic traps (like a quartic plus quadratic trapping potential) and traps with a fixed boundary: 
In the latter case the transition takes place in a parameter regime where the size of vortices is very 
small relative to the width of the annulus whereas in 'soft' traps the vortex lattice persists until the 
width of the annulus becomes comparable to the vortex cores. Moreover, the density profile in the 
annulus where the bulk is concentrated is, in the 'soft' case, approximately gaussian with long tails 
and not of the Thomas- Fermi type like in a trap with a fixed boundary. 

PACS numbers: 03.75.Hh, 47.32.-y, 47.37.+q 



I. INTRODUCTION 



A superfiuid confined in a rotating trap undergoes sev- 
eral phase transitions as the rotational speed is increased. 
In anharmonic traps, where the speed can in principle be 
arbitrarily large, these transitions can essentially be as- 
sociated with three critical speeds. At slow rotation the 
fluid is vortex free [l|, [2[ but when the speed exceeds a 
first critical value a quantized vortex is created. As the 
speed goes further up the number of vortices increases 
[3|48| and a vortex lattice emerges [43J. This picture holds 
in harmonic and anharmonic traps alike, but in the latter 
case a new feature comes into play when a second crit- 
ical velocity is exceeded: The centrifugal forces create a 
'hole' with strongly depleted density around the center 
of the trap [914 llj while the vortex lattice still prevails in 
the bulk. At a third critical speed a remarkable transi- 
tion takes place: The vortex lattice disappears and the 
fluid becomes again vortex free in the bulk. Thus the 
process described above is in a sense reversed. All vor- 
ticity is now concentrated in a giant vortex situated in 
the hole and creating a macroscopic phase circulation in 



the bulk. In the past several authors have studied this 
phenomenon theoretically by variational and numerical 
methods [9], Il34l8j but mathematically rigorous proofs 
of the giant vortex transition have been obtained only 
very recently [19I423J . An experimental realization of this 
transition appears to be still out of reach although an- 
harmonic traps have been available already for some time 

2U22I. 

In this paper we present rigorous results on the giant 
vortex transition in a two-dimensional trapping poten- 
tial that is the sum of a quadratic and a homogeneous 
potential of the form 



y trap ( r ) = kr s + \nl s y 



(1) 



with r the radial variable, k > 0, s > 2 and < fl osc < 
£! ro t where r2 rot denotes the rotational speed. The case 
s = 4 was studied in [19| in an asymptotic regime that 
corresponds to a fixed value of the interaction strength 
while the rotational speed tends to infinity. In contrast, 
the papers 20H23J focus on the interplay between rota- 
tion speed and interaction strength and provide precise 
information about the third critical speed as a function 



of the interaction parameter when the latter is large. The 
model studied in these papers is that of a 'flat' trap with 
the unit circle as boundary as in lfJ,LL! | that can formally 
be regarded as the limiting case s — oo. A mathematical 
advantage of this model is that the extension of the sys- 
tem is fixed, while for finite s the system expands as f2 
and/or e _1 tend to oo. There are, however, both physi- 
cal and mathematical reasons for treating the latter case 
separately. One reason is that s — 4 corresponds to the 
lowest correction beyond quadratic in the Taylor expan- 
sion of a symmetric potential around its minimum and 
such a deviation from a quadratic potential has a better 
chance to be realized in experiments than the limiting 
case s — > oo. A further reason is that the limit s — > oo 
can not be interchanged with the limit of strong cou- 
pling which is the basis of almost all rigorous analysis of 
quantized vortices including the present one. This will 
be explained further below. In particular, the formulas 
for the limiting velocities obtained in this paper do not 
simply pass in the s — > oo limit over to those for a flat 
trap with Dirichlet boundary conditions as considered in 
21], contrary to what might be expected. 

In fact, our analysis reveals significant qualitative dif- 
ferences between the two cases. In the flat trap the gi- 
ant vortex transition takes place at a rotational velocity 
where the vortex cores are still vanishingly small rela- 
tive to the width of the annulus containing the bulk of 
the density. In a trap with finite s, on the other hand, 
the vortex lattice persists until the width of the annu- 
lus becomes comparable to the size of the vortex cores. 
Moreover, the density profile in the flat trap is well ap- 
proximated by a 'Thomas-Fermi' (TF) functional with- 
out a kinetic energy term, while in the case of finite s 
the radial kinetic energy can not be neglected and the 
profile in the radial variable is approximately gaussian 
centered at the middle of the annulus. The long tails of 
such a function impose the use of a larger domain than 
in former situations 20j, |2lj, [28 1 . All together these dif- 
ferences necessitate new ideas for the proofs of the giant 
vortex transition compared to the earlier papers. Com- 
mon to the present setting and 20|, |21| is a macroscopic 



phase circulation around the annulus as well as breaking 
of rotational symmetry of the density in the ground state 
even in the giant vortex regime. 

The mathematical proofs of some of the statements in 
the sequel are rather lengthy and will not be detailed in 
the present paper that is concerned with essential ideas 
and the main results. A full account of the proofs can be 
found in [22j . 



II. THE MATHEMATICAL SETTING 

We now define precisely the mathematical setting 
which is that of two-dimensional Gross-Pitaevskii (GP) 



theory, cf . [3| . The general form of the energy functional 
for the superfluid order parameter "J/ (wave function of 
the condensate) in a rotating trap is 



£%M 



dr 



{W 



iV + A)*| 



+ (^ r ap-itt r V 2 )l*| 2 



1*1 



(2) 



Here A = £l rot e3 A r with 51 rot > the rotation veloc- 
ity, e3 the unit vector in the ^-direction, r = [xi,X2), 
r = |r|, Vtrap the trap potential and 1/e 2 with e > the 
GP coupling constant. The latter is, for a dilute Bose 
gas, given by 2irNa/L with N the particle number, a the 
scattering length of the interaction potential between the 
particles and L a characteristic length in the ^-direction 
[29 ] . Units have been chosen such that H and the particle 
mass are both 1. (This differs from [20 . |21 | where the 
particle mass is taken to be ^)- The normalization of the 
wave function is J |\I>| 2 = 1. The subscript "phys" indi- 
cates that the functional ([2]) is written in the terms of the 
original physical variables, in contrast to the scaled func- 
tional defined in © below. We denote by E^ ys the GP 
energy, i.e., the minimum of @ under the normalization 
condition, and by v[/ GP any of the (in general non-unique) 
minimizers. In this paper we always assume strong cou- 
pling, which means e<l. 

We now specialize to external potentials of the form 
([1]) and 51 rot > fi osc > 0. With the definition $7 c ff = 
(fi 2 ot — f2 2 sc ) 1//2 the sum of the external and centrifugal 
potentials in @ becomes 



Hrap(r-) - 5^ 2 ot 7 



kr 



2^cff 



2 2 

r . 



(3) 



The limiting case s — > oo and fL s ^ = corresponds to 
the 'flat' trap considered in lfj, [ll|, [2CJ, [21| ; the potential 



([3]) is then simply — i£7 2 ot r 2 and the integration is limited 
to the unit disc in R 2 . 

In order that the effect of the quadratic term in the 
potential is visible also when J7 rot — > oo it is natural to 
keep the ratio f2 c ff/^rot fixed and we write accordingly 



"off - 7^? 



(4) 



with < 7 < 1 fixed. The potential (kr s — \^Vi^ ot r 2 ) has 
a unique minimum at r — R m with 



Rm — 



V sk 



l/( s -2) 



(5) 



We now write 

r = i? m x, r = R m x, *(r) = i? m V(x), ^rot = -R m 2 ^ 
and obtain 



(G) 



£^ ys m = R m 2 £ GP m (7) 



with the scaled energy functional 

£ GP M= / {iKiv + nze^i 2 

JR 2 

+- f n 2 V(x)\tP\ 2 +e-' 2 \ij\' l }d 2 x (8) 
where we have written 



V{x) = (±x s 



¥ 2 ) 



(9) 



Note that the scaled potential r yil 2 V(x) has a unique min- 
imum at x = 1, independently of f2, while the minimum 
of ([3]) wanders to infinity as J7 rot — > cxd. Note also that 
we can take Q — > oo either by letting J7 rot — > oo keeping 
fc fixed, or by taking k — > at fixed ri rot > 51 osc . We 
note further that by ([5]) and © the original rotational 
velocity J7 rot is related to $7 by 



rU = (*fc/7) 2/(s+2) ^ 2)/(s+2) - 



(10) 



In particular, for the important special case s = 4, 

n rot - fi 1 / 3 . (ii) 

The potential term and the interaction term in (|5J) be- 
come comparable when £1 ~ e _1 . As discussed below, 
this is the order of the second critical speed Q C 2 above 
which the centrifugal force creates a hole. We are pri- 
marily interested in the case of fast rotation well above 
the second critical speed which means that f2 » e _1 . 

For tt < e _1 it is more convenient to use a dif- 
ferent scaling tha n (EH , replacing R m ~ ^ r ot by 
Re ~ e- 2 /( s + 2 ) [33, Ej. In terms of the scaled rotational 



velocity f2' = R 2 Vt rot the first critical velocity, where vor- 
tices start to appear, is Q,' cl ~ |lne|. See [IIlZI, [8( for the 
case of harmonic traps and [31] for an adaption to s > 2. 
Note that ft' ~ £ _1 is equivalent to fl ~ e^ 1 . 



£~ 4 the radial kinetic energy significantly influences the 
bulk density profile and the TF approximation becomes 
inaccurate. This will be discussed further below. 

From (fT2|) it is clear that p TF vanishes at the origin 
for p TF = and a hole of fine radius forms as soon as 
/i TF < 0. The normalization of (fT2"|) implies that the 
critical velocity for the appearance of the hole is given 
by 



n c2 = s- 1 ((2/ 7 )/[-V] + ) 



-1/2 



(14) 



As (eft) ->• co we have ^ TF /(7fi 2 ) ->■ (s-2)/2s and the 
density p TF becomes concentrated around a; = 1 . The 
inner and outer radii, a;i n < 1 and x out > 1 respectively, 
of the support, as well as the chemical potential ^t TF , are 
determined by the equations 



p TF (zi„) = p TF (a;out) = 0, 2tt 



p (x) xdx = 1. 



. (15) 
A Taylor expansion of V around its minimum (maximum 

for p TF (x)) at x = 1 (see Section 2.3 in [30J for details) 

gives the thickness of the support: 

xout - x in = ( £ n)- 2 / 3 (12/( s - 2) 7 ) 1 / 3 (1 + 0((ef>r 2/3 ))- 

(16) 
By the normalization of p it follows that the maxi- 



mum ||/3 TF ||oo 



p TF (l) is 0((el2) 2 / 3 ). In the flat trap, 



on the other hand, that corresponds formally to s = 00, 
the thickness of the annulus where the TF density is con- 
centrated is 0((efi) _1 ) and density of order O(efl) [111 ]. 
The reason for the different powers of ef2 can be un- 
derstood by the following consideration. The Taylor ex- 
pansion leading to (fl6|) is justified as soon as the turning 
point return where V^"(a;turn) = is much farther from 1 
than the inner and outer radii x- m and x ou t, which means 
that 



III. THE TF DENSITY PROFILE 

In the parameter range Q <C £~ 4 the bulk density pro- 
file of a minimizer ip of ((5J) can be approximately de- 
scribed by the Thomas-Fermi (TF) density 



P TF (*) = y[M TF -7tt 2 ^ 



(12) 



where [t]+ = t if t > and zero otherwise. The chem- 
ical potential p TF is determined by the normalization 
f p TF = 1. The density p TF is the minimizer of the TF 
functional 



£ TF [p] = / [ 1 Vl 2 V(x)p{x) + e- 2 p(x) 2 } d 2 x, (13) 

JR 2 

i.e, the GP functional ((SJ) without the kinetic term. The 
corresponding energy will be denoted E TF . For fl > 



l-a;tur„»(£tt)~ 2/3 (s-l 



-1/3 



(17) 



Now xtum = (l/( s — l)) 1 '^ ^ and since we are in- 
terested in large s we can write (JTTJ) as 1 — s _1 / s 3> 
(eft)-2/3 s -i/3. Since 1 _ s -i/s = 1 _ ex p(-(ln s)/s) = 

(In s)/s + 0((lns/s) 2 ) we obtain 



efi>s/(lns) 3 / 2 



(18) 



as condition for the validity of the Taylor expansion. 
While this condition is always fulfilled for each finite 
s if efl is large enough it is clearly violated for every 
fixed value of efl if s — > co. In fact, in the flat trap the 
Taylor expansion fails and the TF density has the form 
pZ(x) ~ (rfJ)V ~ x 2 n } + (cf. 11], Eq. (A.7)) with 
maximum value p^ F t (l) ~ (efi). 

In the following we consider a fixed, finite s and employ 
formula (fH)|) above. 



IV. THE VORTEX LATTICE REGIME 

An upper bound for the ground state energy E GP of 
(J8]) can be obtained by a variational ansatz that is analo- 
gous to Eq. (4.1) in [llj. It corresponds to a bulk profile 
determined by the TF density and a regular lattice of vor- 
tices localized at positions Xj in the disc with radius x ou t 
centered at the origin. More precisely, the trial function 
is of the form 



•*K x ) = c vpMC( x )0( x ) 



(19) 



where c is a normalization constant, p a suitable regular- 
ization of the TF density p TF , £(x) a function vanishing 
at the lattice points Xj and c£(x) = Ylj(C — CjVIC — 01 
a phase factor generated by vortices of unit strength in 
each lattice point. We have here used the complex nota- 
tion £ = xi +ix2 for points in R 2 . The vortices are placed 
so that V</> compensates as far as possible the vector po- 
tential term proportional to ft in the kinetic energy which 
means an arrangement in a triangular lattice with den- 
sity ft/71- . Moreover, if t is the radius of a vortex core 
where the function £ deviates significantly from 1 the ki- 
netic energy of a vortex localized in Xj is, to lowest order 
in the small parameters, ~ p(xj)| ln(t ft)]. Optimizing 
t to minimize the sum of kinetic and interaction energy 
gives t ~ £•/ p(xj) 1 / 2 , provided t is much smaller than 
the distance between vortices which is ~ ft" 1 / 2 . With 
p(x.j) ~ (eft) 2 / 3 this leads to t ~ e^ft" 1 / 3 and this is 
•C ft -1 / 2 if ft -C £~ 4 . Following closely the computation 
in Section 4 in [ll| one now obtains for e^ 1 < ft <C e~ A 
the upper bound 

0((e 4 ft) 1/3 ). (20) 



E GP < E TF + m\\n(e 4 n)\(l 



The last term is the radial kinetic energy of the density 
profile p. It is smaller than the second term if ft <C £~ 4 . 
A lower bound matching (|20|) is considerably more dif- 
ficult to achieve, but it can be proved using techniques 
from Ginzburg-Landau (GL) theory in the same way as 
in Section 5 in [111]. The result is 



Theorem 1 (Energy between ft C 2 and ft C 3) 
Ife^ 1 < ft < e~ 4 as e -)■ ; then 



E 



GP 



E 



TF 



±ft|ln( £ 4 ft)|(l + (l)). 



(21) 



An important difference to the flat trap considered in 
11 1 becomes apparent here: In |11| an upper bound cor- 
responding to (I2TJ1) (with ln(e) in place of ln(e 4 ft)) is 
shown to be valid under the condition ft <C £~ 2 but the 
lower bound, derived by using techniques from GL theory 
[H, holds only for ft < £- 2 | lnel" 1 . 

For rotational speeds between the first and the second 
value the energy asymptotics can be proved in a similar 
way and is given by 



Theorem 2 (Energy between ft c i and ft C 2) 
// | lne| < Q.' < e _1 as e -)■ 0, then 

E GP ' =E TF ' + in'|ln(e 2 0')l(l + o(l))- (22) 

Here ft' = i? 2 ft rot as before while E GP ' and E TF ' de- 
note the GP energy and the TF energy respectively, mul- 
tiplied by R~ 2 rather than i?~ 2 . 

A further result that holds in the regimes of both The- 
orem 1 and Theorem 2 is that the vorticity is uniformly 
distributed in the bulk in the limit e — > 0. A proof for the 
case of a flat trap with Dirichlet (or Neumann) boundary 
conditions is given in 2l|, Theorem 1.3, and can be gen- 
eralized to the present situation. The precise formulation 
of the statement in the regime of Theorem 1 is delicate 
because of the concentration of the density in an annulus 
that gets thinner as e — > 0, but the main point is that the 
phase circulation around a subset S with area \S\ in the 
annulus is asymptotically equal to 2ft |<S|, and this holds 
uniformly in S provided |<S| is not too small. 



V. THE GIANT VORTEX REGIME 

The first step in a study of the giant vortex transition 
is to consider a variational ansatz for the wave function 
of the form 



-0(x) = g(x) exp(ift^) 



(23) 



with a real valued function <?, normalized such that 
j g 2 = 1. The ansatz (|23|) is well behaved as a function 
of the angular variable $ if ft is an integer, otherwise 
it should be replaced by the integer part [ft]. In order 
to simplify the notation, however, we shall in the sequel 
always assume that ft is an integer; since ft — > oo the 
inclusion of the difference ft — [ft] for non-integer values 
leads only to negligible corrections. Inserting (|23|) into 
© gives 

£ GP m= I {iiv^ + ift^-x-w 

Jr 2 
+ 7 ft 2 (ix s - \x 2 ) g 2 + e-y } d 2 x = £^[g]. (24) 

The unique positive minimizer g gv of the functional 
£ gv [g] is rotationally symmetric, i.e., a function of the 
radial variable x alone. The corresponding energy will 
be denoted E sv . A rough upper bound for it can 
be obtained by taking for g a regularizion of yp^F. 
Since g is concentrated in an annulus of width £ — 
0((eft)- 2 / 3 ) and g 2 = 0((eft) 2 / 3 ) the angular con- 
tribution to the kinetic energy, ift 2 J(x 



-1^2„2 



)V, is 

0((m) 2 ) = <3(£- 4 / 3 ft 2 / 3 ) while the radial kinetic term 
/ |V ff | 2 is 0((eft) 4 / 3 | ln(e 4 ft)|) as in flUD. Hence 

£ GP <^ TF +0(e- 4 / 3 ft 2 / 3 )+0((eft) 4 / 3 ]ln( £ 4 ft)|). (25) 



with 



From now on we shall always assume that 

Q = ft a e- 4 (26) 



with some fixed Slo > while e — >• 0. (For the physical 
rotational velocity (|TU1) and s = 4 this means that O rot ~ 
£ -4/3 ) Thcn the seconc i t e rm m (g5]) is 0(O 2/3 /e 4 ) while 
the second term in (j2Tj> (vortex lattice kinetic energy) 
is O(0o| In Slo|/e 4 ) and thus larger if fio is sufficiently 
large. A radial kinetic energy term O(fl ' |lnOo|/e 4 ) 



is common to both (121|) and ([25]) . The bottom line is 
that for large Qq the giant vortex ansatz is energetically 
favorable to (|2Tj) . 

These simple considerations are, however, far from a 
proof that a true minimize!' ^ GP of f5J has no vortices in 
the bulk above some fi c3 ~ e~ 4 . In [2Gj, |21| such a proof 
is carried out in full detail for the case of a flat trap 
(s = oo), both with Neumann and Dirichlet boundary 
conditions, and it is shown that there C 3 ~ e~ 2 | lne| _1 . 
The technique used in that proof depends on tools that 
were originally developed in the context of GL theory, 
in particular vortex ball constructions and jacobian es- 
timates [32H35J . A prerequisite for these techniques to 
apply is that potential vortices can be isolated in small 
discs with radius much smaller than the thickness of the 
annulus where the bulk of the density is concentrated. 
As pointed out in the discussion preceeding Eq. ([20|) 
above, the radius of vortices is expected to be of order 
e 2 / 3 f2 -1 / 3 ~ £ 2 £V while the thickness of the annulus 
defined by the TF profile is - (eO)- 2 / 3 - e 2 tt~ 2/3 . It is 



thus clear that the methods used in 2C], |21| to prove the 
transition to a giant vortex in a 'flat' trap do not apply 
in the present situation. Nevertheless the absence of vor- 
tices in the bulk can be proved for J7 as in ([26)) provided 
Oo is sufficiently large. The rest of the paper is devoted 
to a precise statement of this result and an outline of its 
proof. 

VI. THE GAUSSIAN DENSITY PROFILE 

In contrast to the regime O <C £~ 4 and the situation 
discussed in 2fJ, [21| the TF profile is not a good approx- 
imation to the bulk density profile in the homogeneous 
trap beyond the vortex lattice regime, i.e., in the present 
situation for $7 ~ e~ 4 . In fact, the bulk of the density is 
contained in an annulus determined by a gaussian density 
distribution that we consider first. We write the energy 
functional (EH) as 



£SV[9] = ~ 7 ^2^ 2 

+ / {||Vfl| 2 + n a tf(aO0 2 +e-y}d 2 * (27) 

Jm 2 



U(x) = I(x-or- 1 ) 2 + 7 & s - \x 2 )+ 1 {s-2)/{2s). (28) 

Taylor expansion of U around x = 1 (e.g. for 1/2 < x < 
3/2) yields 



U{x) = W{x-lf + 0{{x-lf) (29) 



with 



a 2 =4 + 7(5-2). 



(30) 



We consider now for fl as in (f26|) the auxiliary one- 
dimensional functional 

£ aux [/] = / {ll/'l 2 + h&a 2 {x ~ Iff + e- 2 / 4 } Ax 

JR 

= fl / {||/'| 2 + i«V / 2 + Oo 1/2 / 4 } dy (31) 

where the variable transformation y = 1 ' 2 (x — 1), 
f(y) = il _1 / 4 /(a;) has been employed. It is clear that 
all three terms in pip are of the same order of magni- 
tude but the importance of the last term diminishes with 
increasing £7 . Without the last term the minimizer is 
the gaussian 



/o S c(y)^- 1/4 a 1/4 ex P {-W}. 



(32) 



In [22J it is proved that the unique positive minimizer 
5 gv of (|27p is well approximated by the minimizer of (|31[) 
and that the latter is, indeed, approximately gaussian for 
large Slo, so that 

ffgv(x) w g osc (x) = 9}' A f osc {9}l 2 {x - 1)). (33) 

In particular, the integral of g 2 over an annulus 



A, = {x : 1 - 0~ 1/2 77 < x < 1 + £T 1/2 r 1 } 



(34) 



tends to 1 if and only if 77 — » 00, even though Qr x l 2 T] -> 0. 
Furthermore, the same holds for the density l^ 013 ! 2 of the 
minimizer of the full GP functional ([8]) , as shown in [22| . 
Thus any annulus of the form (f34|) contains the bulk of 
the density if r\ — > 00. For the proof of absence of vortices 
in Aiq it is, however, necessary to restrict r\. In fact, we 
prove that the annulus is vortex free if 77 = 0(\ lne) 1 / 2 ). 
In the course of the proof slightly larger annuli, with 
T) = 0(|lne| 3 / 2 ) and 77 = 0(|lne| 3 / 4 ) respectively, have 
also to be considered for technical reasons. 



VII. ENERGY ESTIMATES AND ABSENCE OF 
VORTICES 

Our result on the giant vortex transition is as follows: 



Theorem 3 (Absence of vortices in the bulk) 

There are constants < VLq < oo and c > such that 
for n = il /e 4 with fl > fj and e sufficiently small the 
minimizer ip is free of zeros in the annulus ,4buik = 
{x: |l-<c| ^cfi^VllnsI 1 / 2 }. 

An essential part of the proof is the derivation of the 
precise energy asymptotics in the giant vortex regime: 

Theorem 4 (Energy in the giant vortex regime) 

For il — flo/e with Qq > 0,q the ground state energy is 



E GP = E^ + 0(| lne| 9 / 2 ) = - 7 ^^>fi 2 

2s 



fl 



n 



1 



2?r V 2ttD, 



n 



+ o(ft- 3/4 ) + oar 1 / 2 ) 



0(|lne| 9 / 2 ). (35) 



An upper bound to the energy is obtained by taking 
in (|27|) a trial function built from the gaussian g OS c(x). 
The lower bound is considerably more delicate and is dis- 
cussed further below. As for the comparison with (|2~Tj) we 
note that the negative first term in (|35[) is the same as the 
leading term in E TF , namely the potential energy in the 
minimum of ^VftVix) at x = 1, while the term propor- 
tional to fl is smaller than the term fi| lne 4 fi| = 51 In J7o 
in (|2~Tj) for large Hq. 

For technical reasons we consider besides the func- 
tional (f27|) also a functional £| v defined by the same 
formula except that the integration is restricted to an 
annulus A v with n = 0(|lne| 3 ' 2 ). Its unique positive 
minimizer g^ can be shown to be close to g osc on the 
smaller annulus A 



Vn ■ 



g ri (x) = (l + O(n 1/ ' i ))g osc (x) 



-l/4x 



(36) 



The corresponding energy £& v (g v ) is denoted by E!= v . 
The choice of rj = 0(|lne| 3 / 2 ) is to some extent arbi- 
trary but the method of sub- und supersolutions [36] used 
in the proof of ([36]) and of the exponential smallness of 
\ip GP \ outside of A^ (that is needed for the energy esti- 
mates) requires that | lne| <C r\ <C £ _1 . 

The next step is a decoupling of the energy functional 
that has been used rep eatedly in analogous contexts in 
GL and GP theory 37[ . We define for x e A n a function 
w(x) by writing 



•0^(x) = ff7J (:r)u(x)exp(ifttf). 



(37) 



Since g v is without zeros, the function u contains all pos- 
sible zeros of the minimizer ip GP in the annulus. The 
variational equation for g^ leads to the lower bound 



E GP > E^ 



■ £ vV 



(38) 



with 



£ v [u] = 



f 5 2 {i|V W | 2 -B.J( u )+e- 2 . 9 2 (l-| u | 2 ) 2 } 

J A„ 



where B = fl(x — x 1 )e$ and J(u) = i(uVt 



(39) 
u*S7u). 



The main task is now to estimate the negative term 
involving q 2 B ■J(m). As usual in the context of GP theory 



(see, e.g., 2C],|2l|,|28|) an essential step is an integration 
by parts. Namely, one writes g 2 H = V^F with the dual 
gradient V 4 - = (—d X2 ,d Xl ) and a potential function F. 
In order to employ (|36[) we also restrict the integration to 
A^rj that can be shown to create only negligible errors if 
i] ^> | lne|. If g v \R\ would be exactly symmetrical about 
x = 1 like g osc we could choose F to vanish on both 
boundaries of the annulus A/^j an d integration by parts 
would give 



0JJB.J(u) 



v^7 



FV^Jfu) 



(40) 



v~ 



Moreover, a simple computation employing Q36|) gives 
\F{x)\ < a- 1 (l + 0(f7o 1/4 ) 3 2 (a;), while (V 4 - • J(u)| < 
|Vu| 2 . Thus, because a > 2, the positive first term in 
(j39f integrated over A^ dominates (l40t for f2 large 
enough. 

This reasoning is, however, not rigorous because <7^|B| 
is not perfectly symmetric about x = 1 and if F is chosen 
to vanish on one boundary, e.g., the inner one, it will not 
vanish exactly on the other. The integration by parts 
then creates a boundary term F(R) § x=R 3(u) ■ d£ with 
R the radius of that boundary. To control the circulation 
integral one would like to use a part of positive kinetic 
energy J A g 2 \\7u\ 2 so the first step is to transform the 

boundary integral into two-dimensional integrals. If R = 
R — cfl^ 1 ' 2 with c small and x is a smooth, monotone 
radial function on [R, R] with x(-R) — 0, x(R) = 1 an d 
|V x xl < Of!, 1 / 2 , we can write 



V ± X -J(u)+ hV 1 -^). (41) 



J(w) • d£ -- 

x=R 

The analogous computation holds for R = R + cQr 1 ! 2 
and the interval [R, R] . To bound the integrals in terms of 
J g 2 \ Vu| 2 , however, we need g v to be large in the interval 
[R, R] (or [R, R\). Because g v is, in fact, very small on 
the boundary of the annulus A^ this strategy runs into 
difficulties. 

A way out is to introduce two potential functions, F± 
vanishing on the inner boundary of the annulus and F% , 
vanishing on the outer boundary. The former is applied 
below the radius i? m ax where g v has its maximum, the 
latter above i? max . The integration by parts now creates 
two boundary terms of opposite signs at x — i? max , and 
one has to estimate 



[-Fl(-Rmax) — -FM-Rmax)] 



J(u) • ae. 



(42) 



x=R IL 



Since g v is large in a neigbourhood of _R max Eq. (|4"l"T) 
can be put to good use: Taking R — i? max we have 
9^ > Cfl 1 / 2 on [R,R] by <j36j) , and from gl]), using the 
normalization of <7 2 |u| 2 , we obtain 



/ 3(u)-d£ <C I ^|«||V«| 

+ CQ-V* f g 2 \Vu\ 2 
<C(n-^ 2 + S)f g^Vu^ + CS- 1 (43) 



for any 5 > 0. Here and in the following C denotes 
a finite, positive constant that may differ from line to 
line. The difference \Fi(R max ) — i*2(-R m ax)| is estimated 
separately. It is small because <7tj|B| is approximately 
symmetric about x = 1 , and making use of the variational 
equation for g n it is shown to be at most 0{Slorf/ 2 ). 
Choosing S — CsFIq rj~ 3 / 2 with sufficiently small Cs we 
now obtain for all f2o > ^o sufficiently large the crucial 
bound 



/ g 2 n {\\Vu\ 2 -V-J{u)}>-CSllr? 



C 



5 2 |Vu| 2 . (44) 



s/w 



with C" > 0. Replacing the integration domain by A n 
only produces a negligible correction if 77 3> lne|, and 
E^ v also differs from E sv only by small terms. This com- 
pletes the proof of the lower bound in Theorem 4. 

A further consequence of (|44[) . combined with the vari- 
ational bound E& v < and the exponential smallness of 
f/' GP | 2 outside A^/rj, is the bound 



A 



e' 2 gt(l-\u\ 2 ) 2 <CnW 



(45) 



Vv 



on the interaction term for Qq > ^o- This leads to 
Theorem 3 by the following reasoning. Using the vari- 
ational equation satisfied by g n as well as the Gagliardo- 
Nirenberg inequality in a similar way as in [20], Lemma 
5.1., one obtains the gradient estimate 



|Vu(x)|<C £ - 2 +( c2q / 2 > 
with a as in (1301) for all x such that 



(46) 



(47) 



\l-x\<cn- 1/2 e 2 \\ne\ 1 / 2 . 



Here it has been used that as a consequence of (|36j) we 
have 



9r, 



(x) > Ce- 1+{c2a ' 2) 



(48) 



for x satisfying (|47j) . 

We now claim that as e — > 0, |1 — |u(x)|| < | lne|~ a 
holds for all a > on the annulus defined by (|47p . 
provided c < (2/a) 1 / 2 . The proof is by contradiction: 
Suppose that |1 — |u(x)|| > | In e | a at some x and 
a > 0. Then the gradient estimate (|4"6")) implies that 
|1 - |«(x)|| > |lne|- a /2 on a disk of radius Ce 2 "^ 2 "/ 2 ' 
around x. We thus obtain 



-2„4 



5,(1 



|2\2 



>Cs 



-2+c'q 



In el 



-2« 



(49) 



This is a contradiction to (|45p for c < (2/a) 1 / 2 and 77 = 
lne| 3 / 2 . Thus |l-|u(x)|| < | lnel" 1 holds, implying that 
u, and hence also ip , is free of zeros in the bulk defined 
by 63). 



VIII. CIRCULATION AND SYMMETRY 
BREAKING 

The degree (winding number) of the giant vortex 
ansatz (|23|) is clearly $7. This can also be shown to hold, 
to very good accuracy, for the true minimizer ijj in the 
giant vortex regime, ensuring a macroscopic circulation 
around the central hole where the density is strongly de- 
pleted: 

Theorem 5 (Asymptotics for the degree) 

// fl is given by (|26[) with Qq > Qq and R is any radius 
satisfying R = 1 + 0(Q,~ 1 ' 2 ) then as e —} the degree of 
ijj GP around the circle with radius R is £1 + O(£lo\ lne| 9 ' 4 ) 

Indeed, a simple computation, using (J37]) gives 



degree of -0 = + i(27r)~ 



x=R 



u 1 \u\d m (u\u\ 1 ) 

(50) 

and the second term is easily estimated exploiting Eq. 
(1431) and the bound J g 2 \Vu\ 2 < Cti^rf that follows from 
(gl| together with Ef < 0. 

According to the previous Theorems 3-5 the ansatz 
(|23p gives an excellent approximation to the energy and 
the qualitative properties of a true minimizer ip GP if 17 
is large enough. Nevertheless, while (l23|) is an eigenfunc- 
tion of angular momentum and its modulus therefore ro- 
tationally symmetric, this is not the case for a true min- 
imizer: 

Theorem 6 (Symmetry breaking) 

No minimizer ip GF in the giant vortex regime flo > ^0 
is an eigenfunction of angular momentum. 

The indirect proof is very similar to the proof of a cor- 
responding result in 21|, Theorem 1.6, that in turn is 
inspired by Theorem 2 in [38j, and [39(. One assumes 
that i/; GP (x) = /(x)exp(im?) with a real radial function 
/ and n£Z. As a byproduct of the analysis in the previ- 
ous section it can be shown that n = f2 ( 1 + O (e 4 1 In e | 9 / 4 ) . 



Moreover, / is exponentially small w.r.t. e outside the an- 
nulus A./jj with r; = | lne| 3 ' 2 . It has a unique maximum 
at x m , close to 1. Define 



A(x) = X f {x) for X E A^/rj, x < x m 

and A(x) — for x > x m as well as a smooth interpola- 
tion to zero inside the inner boundary of A/fj. Likewise 
define 



B(x) — nxf(x) for x G A 



s/v 



with a smooth interpolation to outside A/^. Consider 
then 



w(x) = (A(x) + B(x)) exp(i(n + !)■&) 

+ {A(x)-B(x))exp(i(n-l)'&). (51) 

A computation, employing the variational equation for 
/, then shows that the second variation of the GP func- 
tional, i.e., the quadratic form Q given by Eq. (2.3) 
in [381 ] . is negative when evaluated on w, implying that 
f(x) exp(i?r#) can not be a minimizer. 



IX. CONCLUSIONS 

We have analyzed the change in the density and vortex 
patterns of a superfluid in a rotating, anharmonic trap 
as the rotational velocity and the interaction parameter 
both tend to infinity. In particular we have shown rig- 
orously that the fluid undergoes a transition into a giant 
vortex state where there are no vortices in the bulk if the 
rotational velocity exceeds a certain limit depending on 
the interaction strength while a macroscopic circulation 
remains. In the paper we have focused on 'soft' trapping 
potentials, e.g., the sum of a quartic and a quadratic po- 
tential, where the problem turns out to differ markedly, 
both physically and mathematically, from the previously 
considered case of a flat trapping potential with a fixed 
boundary. The differences concern both the shape of the 
bulk density in the giant vortex state, that in soft traps 
turns out to be approximately gaussian rather than of 
Thomas-Fermi type, as well as the relative size of vortex 
cores and the annulus where the bulk is concentrated. 
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